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ABSTRACT 

The abundance of rich clusters is a strong constraint on the mass power 
spectrum. The current constraint can be expressed in the form u^Qn = 0.5 ±0.1 
where ag is the rms mass fluctuation on 8 h^ 1 Mpc scales, fl m is the ratio 
of matter density to the critical density, and 7 is model-dependent. In this 
paper, we determine a general expression for 7 that applies to any models 
with a mixture of cold dark matter plus cosmological constant or quintessence 
(a time-evolving, spatially-inhomogeneous component with negative pressure) 
including dependence on the spectral index n, the Hubble constant h, and 
the equation-of-state of the quintessence component w. The cluster constraint 
is combined with COBE measurements to identify a spectrum of best-fitting 
models. The constraint from the evolution of rich clusters is also discussed. 

Subject headings: cosmology: theory - dark matter - large-scale structure of universe 
galaxies: clusters: general - X-rays: galaxies 
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1. Introduction: 

Rich clusters can be used to constrain cosmological models of large-scale structure 
formation. Rich clusters are the largest virialized objects in the universe and, hence, their 
abundance and evolution can be simply related to the linear mass power spectrum, P(k). 
Their x-ray temperature can be used to infer the cluster mass. Then, Press-Schechter 
(1974) theory can be used to relate the observed cluster abundance as a function of mass to 
the rms fluctuation on 8 h~ x Mpc scales, as- The result is a constraint on a combination of 
parameters, cr&VL^ where 7 is a function of model parameters. 

In this paper, we obtain a general expression for 7 which applies to a wide range of 
models including the standard cold dark matter (sCDM) model, models with a mixture 
of cosmological constant (A) and cold dark matter (ACDM), and QCDM models with 
a mixture of cold dark matter and quintessence, a dynamical, time-evolving, spatially 
inhomogeneous component with negative pressure. An example of quintessence would be 
a scalar field (Q) rolling down a potential V(Q). It has been shown that a spectrum of 
ACDM and QCDM models satisfy all current observational constraints (Wang et al. 1998). 
Our expression for 7 includes dependence on the spectral index n, the Hubble constant h 
(in units of 100 km sec -1 Mpc -1 ), and the equation-of-state of the quintessence component 
w. 

In Section 2, we discuss our derivation of the cluster abundance constraint, which 
follows earlier derivations to some degree but includes some new features when applied 
to quintessence. Some may wish to proceed directly to the resulting constraint, given in 
Section 3. Our result appears to differ slightly from some earlier works, but we explain in 
this section the reasons for those differences. Our purpose in determining 7 is to transform 
observations into a powerful constraint on models. In Section 4, we discuss the further 
constraint derived from studying the evolution of cluster abundance from redshift z = to 
z ~ 1. In Section 5, we conclude by applying the cluster abundance constraint on a 8 in 
combination with the constraint on COBE normalization of P(k) to pick out a spectrum of 
best-fitting ACDM and QCDM models. 

2. Derivation of the Cluster Abundance Constraint for Quintessence Models 

In this section, we present the derivation of the cluster abundance constraint on cr§ for 
QCDM models. We take a pedagogical approach in which we first discuss each step for a 
ACDM model and then point out the differences that arise in QCDM models. The final 
result can be found in Section 3. 
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2.1. Mass-temperature relation 

Press-Schechter theory relates the number density of clusters to their mass. 
Observations determine directly the temperature instead of the mass of clusters. Therefore, 
the mass-temperature relation derived from the virial theorem (Lahav et al. 1991; Lilje 
1992) is used to apply the Press-Schechter relation. 

First, let's consider a universe with vacuum energy density p A . For a spherical 
overdensity, we have 

T mr = ~U G + U A (1) 

where T vir is the kinetic energy, Uq is the potential energy associated with the spherical 
mass overdensity, and U\ is the potential energy associated with A. The kinetic energy 
is T vir = ^MV^ ir where M is the mass of the cluster and V? ir is the mean square velocity 
of particles in the cluster. The gravitational potential is Uq = — f^ff - where G is the 
gravitational constant and R is the radius of the cluster. The potential due to vacuum 
energy is U\ = —^nGp^MR 2 . Re-expressing the background matter energy density at 

redshift z as pb = g^fi (l + z) A where H the is the Hubble parameter, Q is the ratio 
of the matter density to the critical density, and z is the redshift, the virial relation (|1|) 
becomes: 

K4^( G ^ 0V ^) 2/S (l + ,)(l-AM|) (2) 

where A c = Pdusterl Pb is the ratio of the cluster to background density, and £Ia{z) and 
Q m (z) are the density parameters for vacuum energy and matter at redshift z, respectively. 
The observed temperature of the gas is 

where pm p is the mean mass of particles; af os is the line-of-sight velocity dispersion; j3 is 
the ratio of the kinetic energy to the temperature and ks is the Boltzmann constant. The 
mass of the cluster is then: 



M 



knT \ 3/2 
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2 n A ( z ) 

A c Q m (z) 



1 -3/2 



h- l 10 15 M Q (4) 



where fp = f u p/ '/3 where f u is a fudge factor of order unity that allows deviation from the 
simplistic spherical model. Using this relation, the mass-temperature relation of a virialized 
cluster can be computed at any z. However, if we evaluate at redshift z = z c , twice the 
turn-around time (that is, t(z c ) = 2t(z ta ) where z ta is the redshift at which the cluster turns 
around, R = 0), then A c (z c , O m ) = A c (Q m ) becomes a function of Q m only (for A and open 
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universes). Note that z c is the redshift at which the cluster formally collapses to R = 
according to an unperturbed spherical solution. 

In quintessence models, the principal difference is that the energy density in Q 
decreases with time, whereas vacuum energy remains constant. The Q-component does not 
cluster on scales less than 100 Mpc (Caldwell, Dave, & Steinhardt 1998). Consequently, the 
only effect of Q on the abundance of rich clusters with size less than 100 Mpc is through 
its modification of the background evolution. We will restrict ourselves to cases where the 
equation-of-state w is constant or slowly varying. In this case, the ratio Vt\(z) /Vt m (z) above 
can be replaced by (1 + z) 3w Qq(z = 0)/f2 . 

For quintessence models, the ratio of cluster to background density, A c \ Zc = A c (Q m ,w), 
is a function of two variables: 

A / \ Pcluster(Zc) >- ( Rta \ 3 / 1 + Zta \ 3 / c \ 

A c {z = z c ) = - =n——\ — — - . (5) 



Pb{z c ) \R v i r J V 1 + Z, 

where ( = p c iuster(zta) / Pb(zta)', Rta and R V i r are the radius of the cluster at z = z ta and 
at virialization, respectively. (The second equality utilizes the standard assumption that 
the cluster has virialized at z — z c .) The factor ( has been computed by solving for the 
evolution of a spherical overdensity in a cosmological model with constant w (see Appendix 
A). We find that C is weakly model-dependent: C = (37r/4) 2 fi m °- 79+0 - 26 ^- aofki, U a for 
— 1 < w < 0. The fact that this expression is weakly model-dependent means that we can 
also apply it to models with time- varying w. 

By the virial theorem and energy conservation, we have: 

\u G {z c ) + 2U Q (z c ) = U G (z ta ) + U Q (z ta ) (6) 



This leads to an approximate solution 

Ryir 1 ?]v / 2 



R ta 2 + th- 3 Vv /2 



(7) 



where r] t = 2( l VL Q (z ta )/VL [Zta ) and Vv = 2C 1 ((1 + z c )/(l + z ta )fVt Q (z c )/Vt m (z c ). The 
fact that r] v is time-dependent if w ^ —1 is because the energy within the cluster is not 
exactly conserved and the temperature is shifting after virialization due to the change in 
the Q-energy, pg, within the cluster. However, since the matter density p rn is much larger 
than pq in a collapsed cluster, this correction is negligible. 
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2.2. Density Dispersion 

The dispersion of the density field on a given comoving scale R is 

POO rl h 

a 2 (R,z) = / W 2 (kR)A 2 (k,z)^f 
Jo k 

where 



(9) 



and 

A 2 (k,z) = 4nk 3 P 5 (k,z). (10) 

P 5 (k) = \5 k \ 2 is the power spectrum and 5t is the Fourier transform of the fractional density 
perturbation (Sp/p) 2 = (S 2 (x)) 

6 * = 6 ^\) = T^p / d 3 x5(x)e^. (11) 

For constant or slowly varying w < —1/3, the BBKS approximation to the power spectrum 
(Bardeen et al. 1986) is reliable. However, if w > —1/3 or if w is rapidly varying, we find 
no general fitting formula; instead, P(k) must be obtained numerically. 

The rms mass fluctuation <Jg(z) = a(R = 8h~ 1 Mpc, z) can be expressed as 
a$(z) = g(z)as(z = 0), where g(z) is the growth factor. The growth factor is proportional 
to the linear density perturbation S p = Sp/p and normalized to g(z = 0) = 1. We find that 
a good approximation to the the growth index is given by (see Appendix B) 

U = Vr^ = n ™^ a ( 12 ) 

a ma 

where a = 1/(1 + z) is the scale factor and 

125 (l-6w/5) 3 V ; U ;; K J 

The growth factor g(z) is obtained from the integral expression 

nl da 



g(z) ~ a exp 



J a CL 



(14) 



We tested the expression for g obtained from Eq. (|i~3D against the value obtained by 
numerically integrating the density perturbation equations; for 1 — Q m between zero and 
0.8, the accuracy is better than 1%. 
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2.3. Press-Schechter Theory 



According to Press-Shechter theory, the comoving number density of virialized objects 
with mass M = 47i\R 3 pb/3 is: 



dn(M,z) 



Pb 



5 C R da(R,z) 



7T M 2 3a 2 (R, z) dR 



exp 



2a 2 (R,z) 



dM 



(15) 



where 5 C = pu ne ar/ Pb is the perturbaton that would be predicted for a spherical overdensity 
of radius R and mass M according to linear theory. Given the observed number density An 
within a certain temperature range AT, Eqs. (|J) and fll5|) can be used to determine the 
normalization of the mass power spectrum cr 8 . 

The major uncertainties in this method are the observational error in the number 
density An and the systematic error in determining the model parameters 5 C in Eq. (|I5|) 
and fp in Eq. (ffl). Specifically 
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(16) 



where u\ = [5 2 /a(R) 2 — 1] , u 2 = — d In a (R)/d In R and u% = —d\na(R)'/d\nR are 
positive and of order unity in the range of interest. By studying spherical models, we find 
that 5 C varies slowly as a function of Q m , 1.6 < 5 C < 1.686. We also find that fp does not 
depend on the cosmological model and can be determined by numerical simulation. 

According to Eq. (f|), the virial temperature corresponding to a given virial mass 
depends on the redshift at which a cluster is virialized. Therefore, to get the number 
density of clusters of a given temperature range today, we need to find out the virialization 
rate and integrate from z = to z = oo. Assuming that the merger of clusters is negligible, 
the Press-Schechter relation, Eq. (|T5|) can be re-expressed as: 



where x 



dn(T, z) 
dz 

5 c /a(R, z). 




Pb 



d \na(R, z) d Inag(z) 



2ttM(T,z)T dlnR 



dz 



x(x 2 



i)e~dT 



(17) 



Lacey & Cole (1993,1994) and Sasaki (1994) have estimated the corrections due to 
cluster merging (see also, Viana & Liddle 1996). The corrections are small. In our results, 
we average the two estimates of the merging correction to obtain our final result. 
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3. Generalized Cluster Abundance Constraint 

The cluster abundance constraint on <t 8 is obtained by comparing the theoretical 
prediction discussed in the previous section to observations. The observed X-ray cluster 
abundance as a function of temperature was presented by Edge et al (1990) and Henry & 
Arnaud (1991, hereafter HA). After a recent correction (Henry 1997) to HA, the two results 
agree. We have fit the theoretically predicted number density vs. temperature curve (the 
temperature function) to the HA data. 



3.1. Principal Results 

Our results can be fit by 

a 8 = (0.50-0.16)fi- 7(n ' e) (18) 

where G = Q n + Qh — (n s — l) + (h — 0.65) where n s is the spectral index of primordial energy 
density perturbations and h is the present Hubble constant in units of 100km s _1 Mpc -1 . 
For QCDM models with equation-of-state w (including ACDM with w = —1), our fit to 7 is 

7(fi, 9) = 0.21 - 0.22w + 0.33fi + 0.250. (19) 

For open models 

7(0,6) = 0.33 + 0.35^ + 0.206. (20) 



3.2. Comparison to Previous Computations 

Many groups have presented similar constraints on as for ACDM and OCDM models. 
In general, all of them are in reasonable agreement with one another and with our result. 
We identify below the sources of the discrepancies, some real and some only apparent, 
when compared to White, Efstatiou, & Frenk (1993, hereafter WEF); Eke, Cole, & Frenk 
(1996, hereafter ECF); Viana and Liddle (1996, hereafter VL); Pen (1997, hereafter Pen); 
Kitayama and Suto (1997, hereafter KS). 

(1) As we argued, an integration (with merger correction) of Eq. (|17D is necessary since 
the mass-temperature relation is redshift dependent. Most groups only applied Eqs. (f|) 
and (|Uj) at redshift z = to obtain their main results, which leads to an overestimate of 
as by as much as 10%. However, some groups (ECF and Pen) fit the number density and 
temperature relations to numerical simulations to normalize their coefficient fp in Eq. ([|). 
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In so doing, they effectively incorporated the integration correction into the coefficient fp 
for the cases that were numerically tested. Consequently their fitted value of fp does not 
represent precisely f u ^/P as defined in Eq. ([|). We shall call their corresponding coefficient 
fp to emphasize that the physical meaning of this coefficient has been modified to include 
intregration over redshift. Since the contribution of the integration is not proportional to 
fp in all cosmological models, it is more precise and physically meaningful to do separately 
the redshift integration. 

(2) The shape of the theoretical cluster temperature function, Eq. (|T7D, does not agree 
equally well with observations for all parameters. The fit is particularly poor for models 
with large T = Q m h and positive tilt (n s > 1). To handle this problem, some groups 
(WEF; Pen; VL) only fit the observed number density at one particular temperature; this 
introduces some arbitrariness and leads to much larger uncertainties depending on which 
temperature is chosen. In our analysis, we fit the theoretical temperature function to all 
the data points provided by HA. 

(3) The recent correction to the HA data results in a correction to the VL results of 20%. 

(4) Most groups (WEF; ECF; VL; Pen) assume a fixed "shape parameter" T = Q m h. We 
found that no single T is valid for all QCDM models. Instead of expressing our results in 
terms of fixed T, we fix h = 0.65. We include the dependence on h (0/J explicitly. 

(5) Most recent analyses (ECF; VL; Pen; KS) adopted similar modeling of the mass- 
temperature relation. However, there is still about 10% to 20% disagreement on the value 
of fp in Eq. (0j) due to the uncertainties of the numerical simulations. We found that the 
most extensive simulation results, those of ECF and Pen, agree very well with each other. 
By normalizing our theoretical calculations to their simulations, we found fp = 1.1. Notice 
that their reported values for the coefficient are fp = 1.4 and fp = 1.3, respectively. This 
discrepancy is due to rolling into fp the integration effect described under (1). Recent 
simulations by Bryan and Norman (1997) also indicate a similar result, once one corrects 
for their slightly higher value of S c . 

(5) Some groups (VL) used the differential temperature function (the cluster abundance 
within unit temperature interval around a center value T cen ) while others (WEF; ECF; 
Penn; KS) used cumulative temperature function (the cluster abundance with temperature 
above a critical temperature T cri ). These two approaches give similar results because 
the cluster abundance drops exponentially with temperature and the cumulative cluster 
abundance is well approximated by counting the cluster abundance around T cri . We 
compared the results obtained by using the differential temperature function given by HA 
and that obtained by using the cumulative temperature function provided by ECF and 
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found them to be in good agreement. However, the error bar of the latter is much smaller: 
most models were excluded by 95% confidence level by the temperature function fitting. To 
be conservative, we used the former to get our results and errors. 



3.3. Error Estimates 

Eq. ([R]) can be used to estimate the total error for cr 8 . From the present scatter of 
numerical simulation results, fp has about 20% uncertainty, 5 C has about 10% uncertainty 
and another 15% uncertainty comes from the observation. Therefore, the net uncertainty 
quoted in Eq. flT8| ) is about 20% corresponding to 95% confidence level. 



4. Evolution of Abundance 

By applying the same theoretical tools, we can also study the evolution of the cluster 
abundance to obtain further constraints on ag and Q m . The current redshift survey results 
were converted to number densities of clusters with their comoving-1.5 mass (the mass 
within comoving radius R com = l-5h~ l Mpc) greater than a given mass threshold M1.5 
(Carlberg et al. 1997; Bahcall, Fan, & Cen 1997). If the mass profile for the cluster obeys 
M(R) oc RP near R = 1.5h~ l Mpc, and the average virial overdensity is equal to A c as 
calculated in Eq. ([5]), then the virial mass M is related to M L5 by 

M= S 6.97x1^V'mJ '" Ml5 (21) 

Eq. ( flUf) can be used to estimate the number density of observed objects at a given redshift. 
We adopted p = 0.64 as suggested by Carlberg, Yee, & Ellingson (1997). 

The log- abundance as a function of z, log 10 (n(M 15 , z)), is roughly linear as a function 
of z for < z < 1 for the models of interest. A useful parameter to characterize the 
evolution of cluster abundance at redshift < z < 1 is A(Mi.s), defined by 

log 10 (n(Afi. B , z)) « A{M x . b )z + B(M 1A ) (22) 

where n(Mi.s, z) is the number density of clusters with comoving-1.5 mass greater than M1.5 
observed at redshift z. The smaller A{M\.h) is> the stronger the evolution is. By applying 
this analysis to models, we found 

(1) Cluster abundance evolution strongly depends on cr 8 : low er 8 leads to strong evolution. 
This agrees with what Bahcall, Fan, & Cen (1997) have found. This is a general feature of 
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Gaussian-distributed random density peaks. 

(2) Cluster abundance is also sensitive to the equation-of-state of quintessence w: low w 
leads to strong evolution. With the same fl , the growth of density perturbations gets 
suppressed earlier in high w models, therefore, they have a weaker evolution in recent 
epochs (0 < z < 1). 

In Figure 1, we show A(M\$ = 8 x l0 14 h~ 1 Mpc) as a function of a$ for some sample 
models which have been chosen because they all fit current observations well (see discussion 
in following section). We allow a 8 to vary from 0.5 to 1.0 with the COBE normalized a 8 
shown as opaque circles. The current redshift survey data (Bahcall & Fan 1998) indicate 
A(M 1 _ 5 = 8 x 10 14 /i -1 Mpc) ^ -5 at the 3a level (with mean equal to -1.7), which is 
consistent with all six COBE normalized models. 



5. Conclusions 

The cluster abundance and evolution constraints, when combined with future 
measurements of the cosmic microwave background, may be an effective means of 
discriminating quintessence and A models. 

The cosmic microwave background (CMB) anisotropy provides a constraint on the mass 
power spectrum on the horizon scale. For a given model, this constraint from large-scale 
anisotropy as measured by the COBE-DMR satellite (Smoot et al. 1992; Bennett et al. 
1996) can be extrapolated to obtain a limit on a 8 that is completely independent of the 
cluster abundance constraint. In Figure 2, we plot the dependence of a 8 on Qq. For each w, 
a different curve is shown. Along each curve is highlighted the range of a 8 — fig consistent 
with the cluster abundance constraint derived in this paper. Hence, the best-fit models 
are those near the middle of the highlighted regions. These are the same models used as 
examples in Figure 1. 

Near future satellite experiments, such as the NASA Microwave Anisotropy Probe 
(MAP) and the ESA Planck mission, will greatly improve upon COBE by determining 
the temperature anisotropy power spectrum to extremely high precision from large to 
small angular scales. Even a full-sky, cosmic variance limited measurement of the CMB 
anisotropy, though, may not be sufficient to discriminate ACDM from QCDM models. 
There is a degeneracy in parameter space such that, for any given ACDM models, there is a 
continuous family of QCDM models which predicts the same CMB power spectrum (Huey 
et al. 1998; White 1998) It is possible that the data points to a QCDM model, say, which 
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lies outside this degenerate set of models, e.g., a model with rapidly varying w. However, 
if the data points to the degenerate set of models, then it is critically important to find a 
method of discriminating models further. Not only does degeneracy mean that A cannot be 
distinguished from quintessence, but also that large uncertainty in Q m and h. Here we wish 
to illustrate how cluster abundance may play an important role. 

An example of a "degeneracy curve" is shown in Figure 3. Given a value of h for 
any one point along the curve, values of h can be chosen for other points along the curve 
such that the models are all indistinguishable from CMB measurements. The near-future 
satellites are capable of limiting parameter space to a single degeneracy curve. The CMB 
anisotropy also narrowly constrains n s , VL b h 2 , and fl m h 2 . However, even when combined 
with the cluster abundance constraint and other cosmological constraints from the age, 
Hubble constant, baryon fraction, Lyman-a opacity, deceleration parameter and mass power 
spectrum, a substantial degeneracy can remain. Figure 3 includes a shaded region which 
exemplifies the range which these models might allow, based on current measurements 
(Wang et al. 1998; Huey et al. 1998). Because of the uncertainty in cluster abundance 
at z = and other cosmic parameters, the overlap between the degeneracy curve and the 
shaded region allows a wide range of Q m , h and w. 

Cluster evolution offers a promising approach for breaking the degeneracy. Figure 4 
illustrates the variation of A(M L5 ) as a function of w for models along the degeneracy curve 
and inside the shaded region of Figure 3. The variation in A(M 1 £) is nearly 2, corresponding 
to nearly two orders of magnitude variation in abundance at redshift z = 0.5. The range of 
A(Mi. 5 ) is between -3.5 and -5.5 in this case, but this could be shifted upward or downward 
by adjusting cosmic parameters. The point is that models which are degenerate in terms 
of CMB anisotropy are spread out in A(M 15 ). If the measurements can be refined so that 
A(M L5 ) is know to better than ±0.5, then cluster evolution may play an important role in 
discriminating between quintessence and vacuum density and, thereby, determining Vt rn and 
h. 

We have benefitted greatly from many discussions with N. Bahcall and R. Caldwell. 
We also thank P. Bode, G. Bryan, A. Liddle, Y. Suto and P. Viana for useful suggestions. 
This research was supported by the Department of Energy at Penn, DE-FG02-95ER40893. 
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A. Spherical Model 



We study a spherical overdensity with uniform matter density p c iuster and radius R in 
a background that satisfies the Friedmann equation: 



8ttG 



{pb + Pq) 



(Al) 



where a is the scale factor, p^ is the background matter energy density and pQ is the energy 
density in Q. Quintessence doesn't cluster at the interesting scales; the energy density in 
Q remains the same both inside and outside the overdensity patch. Because the curvature 
is not a constant inside the overdensity patch, we use the time-time component of the 
Einstein equations (which does not involve the curvature term) to solve for the growth of 
the overdensity patch 
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We have used pq = wpQ to obtain the second equality. Now, we define 
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where a ta and R ta are the scale factor and the radius at turn-around time, then 
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where fl m (x) is the matter energy density parameter at x, H ta is the Hubble constant at 



turn-around time and ( = (pciuster/ Pb)\x=i- Eqs. (|A1|) and ( |A2| ) can be then written as 

(att m (x))- 1 ' 2 (A* 



dx 
dr 
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where dr = H ta ^jVt m {x = l)dt. With the boundary condition dy/dr\ x =i = and y\ x =o = 0, 
( is uniquely determined by Eqs. ( |A8| ) and ( |A9| ) , given the function form of Q m (x). For 
constant w, we have 

+ 1 a""'^ 1 V T (aio) 

and £ obtained from Eqs. (|A"g| ) and (|A9|) can be well fitted by 

C = f^lY n -0.79+0.26n ra -0.06«,| a|=i (An) 

The linear overdensity 5 C at £(z c ) = 2t(z ta ) can also be calculated by evolving Eqs. ( |A8| ) and 
(|A9|). At early time, the perturbation is linear 

Pduste L= M (A12) 

Pb \yj 

Once 5 C is known at some xq, then it is easily obtained at an arbitrary time 

5 c (x) = 4^U(x ) (A13) 
where g is the growth factor that can be calculated by using Eq. fllij) . 



B. Growth Index in a Quintessense Universe 

The Q-component does not participate directly in cluster formation, but it alters the 
background cosmic evolution. The linear perturbation equation can be written as: 

5 + 2-5 = 4tcG P 5 (Bl) 
a 

where a is the scale factor of the Universe, dot means derivative with respect to physical 
time t, 5 = 5p m /p m , p m and 5p m are the density and overdensity of the matter respectively. 
The background evolution equations in a flat universe are: 

d\ 2 8nG , , . , 

"J =—(pm + P Q ) (B2) 

a ( d s2 



2- + ( - ) = -8ttGw Pq (B3) 

where pq is the energy density of the Q-component and w = Pq/pq is the equation-of-state 
of the Q-component. Now, we can define a matter energy density parameter Q = Q(a) so 
that: 

H 2 = (a\ 2 = 8nGprn m 

\aj 3 
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From Eqs. 



3|) and ( |B4|) , we can get: 

1 3 



aa 
T2" 



w(l-ft) 



By using Eqs. 



By using Eqs. (|BT 



a- 2 2 
and conservation of stress energy d(pa 3 ) 

dm a 
and ( p5|) we get: 



-pd(a 3 



d 2 ln5 /dlnef 



diner 



+ 



din 5 



din a 



The growth index / is defined as: 



din a 
din 5 



1 3 n 
w(l 

2 2 v 



ft) 



5 n 

2 



din a 



By using Eqs. (|B6|) and (|B7|) , we are able to get the equation for / in terms of ft: 



3™ft(l-ft)^ + J 
ail 



1 3 n 
w(l 

2 2 v 



ft) 



+ f 2 = -ft 
17 2 



Now, we introduce variable a, so that / = ft a , and Eq. ( |B9| ) becomes: 

3w(l - ft)ftlnft^ - 3w(a - ^)ft + ft Q - -ft 1 "" + 3wa - |w + i = 

For slowly varying equation-of-state (|dif/dft| <C 1/(1 — ft)), we shall get: 

3 3 (1 -w)(l -3w/2) 



a 



5- 



+ 



1 — 10 



125 (l-6w/5) 3 



;i-ft) + o((i-ft) J ) 



(B5) 



(B6) 



(B7) 



;bs) 



(B9) 



(BIO) 



(bii; 



By following a similar derivation, we found that Eq. ( |B11|) is valid for an open universe if 
we set w = —1/3. Hence, a is weakly dependent on ft m . The result is a ~ 6/11 for ACDM 
(w=-l) and a « 4/7 for OCDM. 
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Fig. 1. — The evolution index of cluster abundance for six models selected because they are 
the best-fit to the combination of COBE and cluster abundance (at z — 0) constraints. The 
model parameters are (from bottom to top): (1) w — — 1, tt m = 0.35, n s — 1, h — 0.65, 
tt b = 0.047; (2) w = -5/6, tt m = 0.34, n s = 1, h = 0.66, tt b = 0.046; (3) iu = -2/3, 
fi m = 0.35, n s = 1, h = 0.66, fi 6 = 0.046; (4) w = -1/2, fi m = 0.36, n s = 1, h = 0.68, 
= 0.043; (5) w = -1/3, Vt m = 0.44, n s = 1, h = 0.67, fi fe = 0.045; (6) w = -1/6, 
tt m = 0.49, n s = 1.1, h = 0.70, tt h = 0.042. 
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Fig. 2. — COBE- normalized erg as a function of £Iq for six constant w models with n s = 1 
and h = 0.65. They are (from top to bottom): (1) w = — 1; (2) w = —2/3; (3) w = —1/2; 
(4) w = —1/3; (5) w = —1/6; (6) w — 0. The highlighted regions indicate where the x-ray 



cluster abundance constraints imposed by Eq. ( |i8|) overlap the COBE constraint. Best-fit 
models correspond to the overlap region. 
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Fig. 3. — Cosmic microwave background (CMB) anisotropy measurements, combined with 
other constraints on cosmological parameters, may be unable to distinguish among a family of 
A and quintessence models, as illustrated here. All models along the degeneracy curve shown 
in the figure produce a temperature anisotropy power spectrum that is indistinguishable, 
given cosmic variance uncertainty. In addition, the CMB can determine other parameters: for 
this illustration, we have assumed fl m h 2 = 1.5, Qf,h 2 = 0.02 and n s — 1 (reasonable values). 
Even if these parameters are determined precisely and combined with other observational 
constraints, there remains substantial uncertainty (shaded region) that may not do much to 
discriminate among the degenerate models, as illustrated here. 




Fig. 4. — The evolution of cluster abundance may break down the degeneracy between A and 
quintessence models illustrated in Fig. 3. v4(M L5 ), the slope of log-abundance vs. redshift 
z, for the models along the degeneracy curve shown in Figure 3. 



